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Abstract. 

Pauli blocking of spontaneous emission by a single excited-state atom has been predicted to be dramatic 
at low temperature when the Fermi energy Bp exceeds the recoil energy The photon scattering rate of 
a ground-state Fermi gas can also be suppressed by occupation of the final states accessible to a recoiling 
atom, however suppression is diminished by scattering events near the Fermi edge. We analyze two new 
approaches to improve the visibility of Pauli blocking in a trapped Fermi gas. Focusing the incident light 
to excite preferentially the high-density region of the cloud can increase the blocking signature by 14%, and 
is most effective at intermediate temperature. Spontaneous Raman scattering between imbalanced internal 
states can be strongly suppressed at low temperature, and is completely blocked for a final-state _Ep > AE^ 
in the high imbalance limit. 
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1. Introduction 



The Purcell Effect is the enhancement or reduction of scattering due to a modification of the electromagnetic 
density of states by a cavity pj. A complementary effect can occur in an ensemble of fermions, where final 
states of the recoiling particle are blocked, reducing the scattering cross section. In semiconductors, such 
blocking creates a Moss-Burnstein shift [5] of the apparent band gap. Further study of this fundamental 
effect has been proposed for neutral Fermi gases [31 131 |5J [51 17] , motivated by the control of trap dimensionality, 
the direct quantification of density, and absence of additional scattering phenomena. 

Light scattering is the primary tool for detection of ultra-cold atoms [8], with a few notable 
exceptions In addition to direct imaging, light scattering has been used to probe density [10], phase 
coherence [TTl [T^], momentum distributions [13 , excitation spectra [TJ, and superradiance [TS] [T5] |T7] 
in quantum degenerate Bose gases. The light scattering properties of degenerate Fermi gases, by 
contrast, have only recently been studied experimentally [T51 |TH], despite numerous theoretical proposals 
[3ll[3[a[2ni[ni[lll[23[Mlll3[Ml[lZl[lHllMllZ]- In-situ optical probes could be particularly useful in 
exploring the physics of paired superfluids [20l [211 [22l [211 [29] . In addition, the temperature dependence of 
light scattering could be exploited for thermometry. 

In this work, we review two scenarios in which Pauli blocking has been considered previously, and then 
consider two new scenarios in which the experimental signature of blocking can be enhanced. Whereas 
pioneering work treated untrapped gases [H [201 1221 [23 > or geometries with spherical [H [23] [23] or cylindrical 
[11 [5] [28] symmetries, we treat a generalized scenario that includes finite temperature and a tri-axial trap 
potential. Since the high optical density of a trapped gas requires off-resonant excitation to avoid the 
multiple-scattering regime, we focus on scattering suppression [H Ul UHl [Z] instead of line shape [31 [SI [?7] . 

We develop, in sections [2] and [3] a semiclassical approach that may be applied to situations in which 
fully quantum calculations have proved onerous. In section [3] we show that our approach reproduces fully 
quantum calculations in the literature. We then calculate angle-averaged finite-temperature signatures 
without imposing any symmetry. Finding that suppression is rarely complete, we evaluate in section [4] 
two approaches to stronger blocking: using a focused excitation beam, and scattering between imbalanced 
populations. Finally, experimental prospects are discussed in section [5] 



2. Methods 



We consider light scattering and spontaneous emission in the presence of a Fermi sea of neutral atoms. 
Atomic excitation is assumed to be far below saturation, and multiple scattering is assumed to be weak. The 
N degenerate fermions are trapped in a three-dimensional harmonic trap with trap frequencies {uji, 102,(^3} 
along three axes. When a subscript is not specified, lu refers to the geometric mean of frequencies. The 
single-particle Hamiltonian is 

^ = J2^p' + l^^'^l (1) 

where m is the mass of the atom, pj is the momentum operator in the jth direction, and qj is the position 
operator in the jth direction. 

We treat spontaneous emission and scattering of an incident photon using a Golden Rule approach, as in 
[H [5] jU The reduction of the scattering rate by Pauli blocking is proportional to the reduction in the number 
of final states, weighted by matrix elements. This is conceptually similar to the suppression of spontaneous 
emission in an optical cavity where the density of electromagnetic states is reduced; here, the density of 
available atomic states is reduced. We define the relative scattering rate S to be the ratio of the scattering 
rate with fermions to the scattering rate with Boltzmann particles, i.e., the rate without blocking effects: 

g,^. ^ E.j:uMv){^-Mu)}\{u\e^''-^\v)\' 
§ This approach docs not treat coherent effects, which are expected to be important within the forward diffraction cone l8l l23l . 
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where k is the recoil momentum of the atom, \v) and \u) are energy eigenstates \vi,V2,V3) and 1*2,1*3), 
and Hi and Uf are the initial and final occupation functions, respectivelj|]|] The matrix element along a 
single direction can be calculated using [32 



V 



2i2 



(3) 



min(uj, Vj), A = 



V 



2 ,.„^2 w\ 

XQj — ^hjlMujj is the ground state width, 
J I, and is the generalized Laguerre polynomial. 



where kj is the projection of fe along the j*^ direction 



Due to the computational demands of a six-dimensional sum, (|2| is more easily calculated in the case 
of a spherically symmetric trap, where three of the sums can be eliminated |3] . A similar approach applied 
to the case of cylindrical symmetry can reduce six sums to five. Without these symmetries, our desktop 
computer did not have the resources to calculate S with experimentally realistic parameters (e.g., N — 10^ 
atoms, and trap frequencies of a;/27r = {500, 800, 30} Hz) at finite temperature. Indeed, no such calculation 
has been published. 

If many states are occupied along all three harmonic axes, a semiclassical integral may capture the 
important physics of the problem. One might also expect a local density approach to be appropriate since 
light scattering predominantly probes high-momentum properties that depend upon local density fluctuations 
in the gas. Starting from the semiclassical phase space element [33] dN = h~^n{t)d^pd?q, where n(e) is the 
quantum statistical occupation function, the relative scattering rate is 



S{k) = 



h 3/ f d^pd^qn,{p,q){l - nf{p+ k,q)} 



(4) 



h^-^ J J d^p d^qni{p,q) 

This local density approach neglects the energy quantization scale set by the level spacing, and thereby 
allows us to rescale the problem into an isotropic one, even though no such symmetry is manifest in the trap 
geometry. The symmetry is broken only by the direction of the momentum recoil due to scatteringj^ A 
semiclassical approximation was also used in [27] to discuss line shape, and to treat the uniform gas in p5]. 
Here we extend these treatments to find the light scattering properties of a trapped gas. 

For large atom numbers and moderate trap anisotropies, we show (in section [3]) that ([4]) is in excellent 
agreement with published calculations based on ([2|. Furthermore, the simplicity of the semiclassical method 

II In Fermi gases, unlike in Bose gases, neglecting the zero point energy when calculating the occupation ni,f leads 

to a fractional error in the chemical potential of (1 + e/2)(6A''e)~-^/'^, where eccentricity e = WzI^jJ^- 

^ For instance, illumination of a cigar-shaped cloud along one of its two radial axes would break all rotational symmetries. 




Figure 1. Comparison of quantum (points) and semiclassical (lines) calculation of spontaneous emission 
blocking. Suppression is absent when S — > 1 and complete when S = 0. For comparison to figure 1 in |4], 
the quantum calculations assume Br = 25?ilo> and spherical symmetry, however these assumptions are not 
necessary for the semiclassical calculation. From top to bottom, Ep/E^i = 0.4, 0.8, 1.2, 1.6, 2.0. The inset 
"cartoon" is a schematic representation of the process under consideration: an atom in the midst of a Fermi 
sea emitting a single photon. 
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allows us to include angular averaging across the scattered photon momentum, to evaluate both finite and 
zero temperatures, and (in section [4]) to consider more complex scenarios. 

3. Signatures of Pauli Blocking 

Two scenarios for Pauli blocking have been considered in previous work. The first is spontaneous emission 
("SE") of a single excited-state atom in the midst of a Fermi sea. The second is light scattering ("LS") off 
a large ground-state ensemble. For reasons discussed further in section 3.1 and in section|5] the LS scenario 



is experimentally more feasible, and is the focus of our work. However, the SE case is the simplest scenario 
that elucidates the blocking effect under discussion. Our treatment can also be compared directly to [4], 
which presents a fully quantum calculation of the same scenario. 

3.1. Spontaneous emission of a single excited atom 

For a single atom in the excited state decaying into a Fermi sea of N atoms, 

ni{e) ^ (Phuj)^ exp{-l3e) and 7i/(e) = [1 + exp/?e]"\ (5) 

where = /sbT, fee is the Boltzmann constant, T is the temperature, and z is the fugacity of the Fermi 
gas. This Ui is normalized by integration over h^^d^pd^q. 

Note that we assume that the excited-state atom sees the same trapping potential as the ground state, 
and is thermalized with the atoms in the ground state. This might be realistic for fermions with long-lived 
metastable states, for instance in rare earth metals j34j in magic wavelength traps [35 . 

Using distributions ([s]), we evaluate Q using a change of variables. The coordinates are rotated such 
that one momentum axis is aligned with the momentum kick fc, and the other five coordinates are combined 
into a five-space radius. We define a dimcnsionlcss momentum ^ = ^/JjEji from the single-photon recoil 
energy Eji = h^k\/2m, where cfc^ is the laser frequency. The normalized scattering rate is then 

SsEitT) = ^z-'e^' dye^^y f^,^{ze-(y+^^'), (6) 

where /„ is the Fermi functiorj^ Note that since we are using a local density method, trap frequencies and 
atom number now affect the scattering rate only through E^. 

+ The Fermi integrals that appear in this section arc of the form 

/■OO 1 

/ '^''-^da—-—-=V{n)U{C), (7) 
Jo e°-/C+l 

where fn{C) is — Li„(— C), and Li„(C) = the polylogarithmic function. 
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Figures [T] and |2] show the essence of the blocking effect: at low temperature, and when the Fermi energy 
is laxger than the recoil energy, the scattering rate decreases. Figure [l] compares our results to Busch et al 
[4] to find that the semiclassical S'se agrees well with a fully quantum calculation, in this case calculated 
using ^ and assuming spherical symmetry. 

Figure |2] shows that complete suppression of spontaneous emission can occur when the Fermi energy 
exceeds the recoil energy. (By contrast, suppression is only complete for infinite Fermi energy when the set 
of initial states is expanded to a full Fermi sea, as discussed in the next section.) At zero temperature, ^ 
takes the simple form 

5sE(«:) = e(K-l), (8) 

where O is the unit step function, and we use a dimensionless momentum k = {Eyi/ E^Y^'^ since ^ is ill- 
defined at zero temperature. The abruptness of this step is the only discrepancy between the quantum and 
semiclassical calculation, as shown at the left-most points in figure [l] Interesting non-isotropic effects have 
been predicted in this limit [4], with possible application to directional single-photon sources Wave- 
function features appear when fceT is comparable to the level spacing hu)j. However, even at the lowest 
observed temperatures of approximately O.OS-Ep/fcB |36ll37j . current experiments remain in the semiclassical 
regime. 

Figure |2]shows that the emission rate 5se('*, T > O.STp) > 0.86, confirming that blocking is not dramatic 
in the non-degenerate regime. However for the weak effects observed at high temperature, one can use a 
series expansion of f5/2' 



5sE(C,r) 



1 - e" 



(-l)"z" 



,e/r. 



11=2 



(9) 



This expression converges for z < 1, i.e., k-QT/E-p > 0.57. 

An approximate form valid for all temperatures can be developed by neglecting the directionality of the 
momentum kick. Given an initial atomic momentum p, the average energy transferred by a kick is simply 
-Er when averaged over a uniform distribution of atomic momenta. Using this energy difference, we can fully 
integrate S: 

SsE{^,T)^z-'ei' Mze-^'). (10) 

Figure [2] compares the approximations ([s]) and (10) to numerical integration of Sse- The average- 
kick approximation ( 10 1 underestimates the blurring of the step function at finite temperatures, but is a 




0.25 0.5 0.75 1 

k(q) 



Figure 3. Scattering rate Sls versus momentum kick K,{a) comparing fully quantum (points) and 
semiclassical (lines) calculations, from Eqs. ^ and | |12| l respectively. Temperatures shown are k^T/Ep 
= 1.0 (triangles), 0.47 (circles), 0.24 (squares), and 0.10 (inverted triangles). The fully quantum calculation 
is for a cylindrically symmetric geometry (lui = lJ2 = '^±, and incident beam along axis of symmetry), and 
parameters are chosen to reproduce figure 3 in [5]: i?R = 21?!Lt;x a^nd aspect ratio is 0.1. Assumptions of 
symmetry, aspect ratio, and Epi/lu± are not used for the semiclassical calculation. 
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Figure 4. Angle-averaged scattering rate Mls versus average kick k for a dipole emission pattern. Thick 
lines are numerical integrals, for k^T/Ep = 0.6, 0.2, and 0.05. The dashed line is the zero temperature 
limit. Inset: Scattering rate Sls versus momentum kick for the same temperatures. 



reasonable estimate at the 20% level and even better at low n. In all cases, suppression is stronger for lower 
recoil momentum (or higher Fermi energies) , since final states fall closer to the centre of the Fermi sea. 

3.2. Light scattering from a large ensemble 

We now consider N polarized fermions in a single ground state, recoiling under the net momentum k of an 
incident and a scattered photon. In the perturbative limit, we ignore the disturbance of removing an atom 
from the distribution, and use the same initial and final state 

n^{e) =nf{e) = [l + z-^exp/3e]~^. (11) 

The initial states have an energetic range that is determined both by temperature and Fermi pressure, unlike 
the case of a single excited-state atom. This makes it easier to scatter out of the Fermi sea, and reduces the 
net blocking effect. 

Integration of ^ using (11) yields a normalized scattering rate 



8 l"^ 1 1 

5Ls(C(a),r) = l--^^y^ a3/2 da j _^ dy ^ ^ ^ ^ ,-i,aHym^)r ' (^2) 

where the results now depend on the angle a between the incident and scattered photon. The rescaled recoil 
momenta are ^(a) = 2^sin^(Q!/2) and K.{a) — 2Ksin^(Q;/2), and we maintain the previously defined quantities 
= [3E^ and = E^/Ep, written without an angle argument. Figure [s] shows a numerical integration 



of ( 12 1 at various temperatures. Points show a reproduction of quantum mechanical calculations assuming 
cylindrical symmetry, for values chosen in [5j. The agreement is excellent. Unlike the SE case, there is no 
disagreement between quantum and scmiclassical calculations at low temperature so long as Ep ^ hujj in 
all directions. 

Since the excitation and decay both contribute a momentum kick, an angle-resolved experiment would 
observe S directly [51 . If an experiment measures the total scattering rate (see discussion in section |5| , we 
observe an angle-averaged suppression factor, which we call M: 



M{k)= S(2kcos^ -)P{a)da, (13) 
Jo 2 

where P{a) = |(1 + cos^ a) for a dipole emission pattern of any polarization, after averaging over the 
azimuthal scattering angle. Numerical evaluation of Mls versus k, no longer a-dependent, is shown in 
figure |4j Comparing Mls and S'ls, shown in the inset, we see that angle averaging produces little qualitative 
change. The low-n limit is identical, but suppression continues to higher k in A^ls- This is due to inclusion 
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of forward-scattering events that produce small kicks and are easy to block. Including these scattering events 
is necessary for quantitative prediction of the scattering suppression. 
As before, we can expand S to find a high temperature expression, 

5Ls(C(a),r) ^ 1 - E E e-«(")^^("-)/", (14) 

however the series converges only for z > 1, where 5ls > 0.91. Keeping only the first term, Si,s{S,{ol)) ~ 
l-f(/3i?F)-Vexp(-e(a)V2). 

At zero temperature, an analytic expression for ^ls can be found: 

SMa)) = 1 - gx(^)e(2 - K{a)), (15) 

X{x) ^ ^^^^^ (-8x4 ^ 26x2 - + ^ co^-\x) . (16) 

A further integral can also be done to find an expression for Mls at zero temperature, and is given in the 
appendix. Both scattering rates are plotted as dashed lines in figure [4] and linearly approach zero as the 
momentum kick goes to zero. 

Quantum corrections to the scattering rate can also be evaluated by considering ^ for various trap 
geometries, and comparing to the geometric insensitivity valid in the semiclassical limit of 3> fujjj for all 
i = {1,2,3} trap frequencies. Evaluating ([2]) for a cylindrically symmetric trap with 20:1 aspect ratio, and 
comparing to a spherically symmetric trap, both for E^ — Whw, we find that Sse changes by less than 5%, 
Sls changes by less than 3%, and M^s changes by less than 2%. 



where 



4. Scenarios with stronger Pauli blocking 

While spontaneous emission can be suppressed fully at finite momentum k for sufficiently low temperature, 
complete blocking is possible only at /c = in the case of light scattering (section |3.2| . The smallest k 
reported to date is « 0.6 |30l [SH |38], so a suppression of at most Mls ~ 0.5 would be expected. In 
this section, we explore two methods to improve the Pauli blocking signal: the use of a focused excitation 
beam ("FB"), and scattering between two imbalanced populations ("IP"). In both cases, we attempt to 
bias scattering toward events with higher local E-p. In the FB scheme, this is done directly by selecting 
the spatial centre of the trap. In the IP scheme, reduced Fermi pressure in the initial state also reduces 




Figure 5. Pauli blocking is enlianced when the excitation beam is focused onto the centre of the cloud. 
The normalized scattering rate Sfb is shown versus beam size a/Rp, for K{a) = 0.5 and k^T = 0.2i?F. 
The arrows indicate the asymptotic limits: the three-dimensional Sls high a, and the two-dimensional 
"-'ls I'l^^t: The same quantities are plotted with a dipole pattern angle-averaged scattering rate 

MpB. 
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Figure 6. The maximum possible strengthening of Pauli blocking with a focused beam, AS = — 5lSi 
is shown versus momentum for temperatures k^T/Ep = 0.5 (cyan), 0.4 (green), 0.3 (red), 0.2 (purple), 
0.1 (blue) and zero (black dashed). Inset: The enhancement AM = M^g — Mls versus normalized recoil 
momentum k, after averaging over a dipole angular distribution. For clarity, fewer temperatures are shown; 
k^T/ Ep = 0.5 (top solid, cyan), 0.2 (bottom solid, purple), 0.1 (middle solid, blue) and zero (black dashed). 
Again the largest effect is at intermediate temperatures. 



the initial kinetic energy. As is shown below, these approaches increase the overlap between accessible final 
states and the Fermi sea of occupation. 

4.1. Focused excitation light 

DeMarco and Jin suggest that stronger blocking might be observed when focusing the incident laser beam 
|5j. Here we evaluate this scheme quantitatively. We consider excitation along a cycling transition, starting 



and ending in the same Fermi sea, as in section 3.2 A focused excitation beam restricts to its intersection 
with the atomic cloud. Assuming the beam propagates along (73, the distribution of initial states is 

91,92) = -^(gi-92)[l + 2:"^exp/3e]"\ (17) 
where I{qi,q2) is a dimensionless intensity distribution of the light. The final state distribution Uf remains 



as before, given by (111 



For simplicity we consider a cylindrically symmetric beam I{qi,q2) = e:xjp {—2ql / a'^ — 2q2/a'^), 
(7 is the waist of the beam. Starting from ^ , we rescale symmetric degrees of freedom and are left with a 
triple integral: 

S'fb (^fa), r, cr) = 1 / da xdx / dy — ttet-, — , , o^ 

^ 1 + 2-le'3'BF(Q+(!/+K(a))2+2:") ' '-^^^ 

where the mean number of atoms excited by the probe is 

N,= ada xdxe-2(-fiW'^)^[i + ^-ie/3£F(a+x^)pi^ (^g) 
Jo Jo 

the radial Fermi radius Rp — ^J2Ey /m/uj_i_, and we have assumed uj±^ = loi — lo-i- Figure [s] shows that 
smaller beam size enhances the suppression. Atoms are excited at the centre of the cloud, where the density 
is higher and thus the local E-p is higher. Since E^i is unchanged, we effectively decrease k. 

In the large cloud (or small beam) limit, a <^ Rp, the spatial selection of the exciting beam becomes a 
delta function. Since rescaled quadratic degrees of freedom are equivalent under the integral, eliminating two 
spatial degrees of freedom is equivalent to eliminating one spatial and one momentum degree of freedom. In 
other words, for a given geometric mean lj, the same scattering rate is observed for a tightly focused beam 
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on an oblate three-dimensional cloud, as would be observed for a two-dimensional cloud with a uniform 
excitation light. This limit is 



da 



dy 



1 



1 



1 _|_ ;^~lf.a+y-^ I _|_ z-^ga+iy+i{a)y 



(20) 



Figure [5] shows as arrows the 3D limit (12 1 and the 2D limit (20 1. 

Figure [b] shows the difference AS = — S'ls (and angle-averaged AM = M^g — Mls) between 
the small- and large-beam limit. This is the maximum effect that changing beam size could have. We 
see that the difference is restricted to AS" < 0.15. Interestingly, the most pronounced effect occurs when 
k-^T / E-p w 0.25. At intermediate temperatures, selecting the centre of the cloud is even more important 
than at zero temperature, since quantum degeneracy varies across the cloud. At lower temperatures and 
momenta, suppression is complete for both the 2D and 3D limits, so focusing is less effective. 

At zero temperature, an analytic expression can be found: 

6(2 -K(a)) 
12^r 



c2D 



K(a)V4- /t(a)2(K(a)^ - 10) + 24cos"i {n{a)/2) 



(21) 



This expression is shown as a dashed line in figure |6] 

A similar expression can be found for the angle-averaged M^g , and is given in the appendix. For a 
variety of temperatures, figures [5] and [6] show the angle-averaged results as insets. In the inset of figure [6] the 
dashed line shows the zero temperature difference between equations (A.l I and (A.2|. As with angle-resolved 
scattering, the enhancement AM is no more than 0.15, and occurs at intermediate temperature. However, 
because of the inclusion of low-a events, suppression is observed (and enhanced) for k > 2. 



4-.2. Imbalanced Fermi gases 

An alternate method of reducing the distribution of initial states is to use the internal structure of the atoms. 
Consider Raman light scattering between two ground states, and the population of atoms split unequally 
between them [36, i37j such that Zi < Zf, where i and / indicate the initial and final states. Now the 
thermalized initial and final distributions are 

n,{e) ^ [1 + z~^ exp I3e]~'^ , and n/(e) = [1 + z"^ exp/3e]"\ (22) 




Figure 7. Scattering rate Sip versus imbalance S, 
temperatures ksT/Ep = 1, 0.5, 0.3, 0.2, and 0.1 
limit from | [24^ . Note that for Figures [t] and [s] the 

quantities, so fc = -Br,/ Ep{. 



at k(o) = 0.5. From top to bottom, solid lines represent 
from | |23[ l; the dashed line shows the zero temperature 
final state Fermi energy has been used for dimensionless 



Enhanced Pauli blocking of light scattering in a trapped Fermi gas 



10 



1 

0.8 
0.6 
0.4 
0.2 









8=0 / 



6=0.5/ 



"f ■ 

1 0.6 
I 

I CQ 
i' 5 0.4 



■ — 



6=0.99' 



0.5 1 

hi 



1.5 2 



0.5 



1 



1.5 



Figure 8. Zero-temperature scattering rate Sjp versus kick K{a), for various imbalances: 5 = 0, 0.5, and 
0.99. This corresponds to balanced, imbalanced, and nearly polarized. Inset: Angle- averaged scattering 
rate Mip versus kick k and at the same three imbalances. Imbalance allows complete suppression to be 
observed for nonzero K. 



As before, we ignore the change in either distribution due to scattered hght or due to interactions. We also 
assume that incident hght excites only atoms from and decays only to |39j . Integration of Q with 
(p2l) yields 



5ip(e(«),T,5) = l 



7r(/3£;Ff)3 7o 



a3/2 da 



dy 



1 



(23) 



where the difference in Fermi energies between the two states is parameterized using 6 ~ (-Epf ~ Epi)/E-p{. 
For instance, Ni ~ Nf{l — S)^, etc. Figures [t] and [s] show the normalized rate of Raman scattering between 
imbalanced Fermi clouds. At finite temperature, figure [7] shows S'lp versus S for k — 0.5, demonstrating that 
an imbalance enhances suppression. Increasing the cloud imbalance enhances blocking because the range 
of initial states is increasingly restricted to lower energies. At low temperature, the effect can be dramatic, 
allowing for complete blocking a,t S > 3/4. 

The zero-temperature limit of Raman scattering is 



16 



16 



Sip{K{a),S) = e{K{a) + b-l){l- —e{b + B)x{-B/h) + —9(1 + b - K{a))x{K{a) + B)/W}, 

OTT 



(24) 



where b = \/l — S, B = {S — k(q;)^)/(2k(q!)), and x is defined in (16). Figure [s] shows zero-temperature 
suppression at various imbalances. In the limit 5 = 0, there is no imbalance and (|24|) becomes (15 1. The 



comparison in figure [8] makes especially clear the wide range of suppression possible at finite recoil momentum 
for imbalanced gases. Complete suppression is possible even in the angle-averaged case for k < 0.5, i.e., for 
> 4£'r. By comparison, in the balanced (or single ground state) case S'ls ~ 0.48 for K{a) = 0.50. 
In the strong polarization limit 5 ^ 1, (24 1 becomes 0(k(q) — 1). This is reminiscent of the spontaneous 
emission case ([8|, apart from the angular dependence of K{a). We can now see new significance in the 



results of section |3.1| in the limit of extreme polarization, the Raman scattering problem is equivalent to 
the spontaneous emission problem with a two-photon recoil momentum. In both cases, suppression is strong 
because a second state allows initial states to be exclusively at the middle of the Fermi sea. 

Comparing the two enhancement scenarios, the FB approach is less effective than the IP scheme, since 
initial state selection occurs only along two coordinates instead of all six coordinatesj^ Both schemes are 
successful at enhancing the expected signal in realistic experimental scenarios: at k = 0.25 and T = O.lTp, 
we find Mls = 0.40, Mfb = 0.32 for a/Rp = 0.50, and Mjp = 0.07 for S = 0.90. The dramatic suppression 

* Another experimentaUy viable situation would be a cigar-shaped cloud with focused excitation in the plane of symmetry. 
This would restrict excitation primarily along qi, for instance, but not q2 or q^, and therefore be less effective than the oblate 
geometry considered in section 14.11 
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that seemed like a distant experimental prospect in section 3.1 (since it required a thermalized excited state 
atom) is feasible using a Raman light scattering scheme. 

5. Experimental realization and conclusion 

To realize a Fermi energy that is I/k^ times the recoil energy, the mean trap strength must be 

For 10^ "^K atoms, oj = 2tt x 860 Hz at Ep = 4Er. For 10^ ^Li atoms, uj = 2tt x 1620 Hz at Ep = 4Er. 

In order to avoid multiple scattering, the atomic sample must have a low optical density. Consider a 
cloud of N fermions at T = in a cylindrically symmetric trap, whose eccentricity is e = w^/wj,. A Thomas 
Fermi profile has a resonant optical density bounded by 

Oi?-- = max{ei/3,e-2/3}3 ( EfV En_ 



Combining with ( 25 1 and the constraint that OD < 1 requires that N be on the order of unity, which is 
clearly incompatible with a Pauli blocking experiment. Quantitative predictions for spontaneous emission 
and for light scattering at resonance will require treatment of multiple scattering. 

Scattering with off-resonant light can avoid this complication. Detuning approximately N^^^ line widths 
away from resonance can achieve OD < 1. Our focus on net scattering rates rather than on line shape is 
partially motivated by this limitation. 

One could use a variety of signatures to search for blocking effects in the lab. Because fewer than one 
photon per atom can be scattered while remaining in the perturbative limit, angle-resolved experiments 
that measure only a fraction of the total scattered light may be difficult. However, a direct measure of 
the integrated rate M, rather than S{a), is possible through absorption imaging. Another measure of M 
would be an optical pumping experiment, in which the efficiency of pumping into an occupied Fermi sea is 
reduced due to blocking effects. In this case the blocking effect would be recorded in atomic populations, 
circumventing lensing effects of the detuned absorption beam by the degenerate cloud. 

In summary, we have presented two light scattering scenarios in which Pauli blocking is strengthened. 
Including the effects of inhomogeneous trapping, finite temperature, and without assuming rotational 
symmetry, we make quantitative predictions for both angle-resolved and angle-averaged signatures. Pauli 
blocking effects can be enahanced by focusing the excitation beam. Dramatic suppression of incoherent 
scattering can be achieved with Raman scattering between internal states for large population imbalances 
and and Ep > 4i?R. Our calculations should aid experimental efforts to observe this fundamental quantum 
optical effect. 
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Appendix: Angle-averaged results for zero temperature 



Substitution of (15) into (13 1 yields an analytic expression for zero-temperature angle-averaged scattering 
in the LS case: 

f 252(5^ {512{VT^ - 1) + 315At (16^2 + 3) sin-^K,) 
Mls(k) = < [2kV1 - ^2(232^4 - 1036^2 + 2157) -t- 2767^1 - - 3456] } for k < 1 (A.l) 

[ 1 - [256 (27k2 + 4) - 9457rK]/(50407rK3) for k > 1 

This quantity is plotted as a dashed line in figure [4] 
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A similar expression can also be found for the FB case in the tightly focused limit: 
2^ { k^ [2V T^ (131 - 32k2) k2 + 263Vr 



336 



64 {Vl 



1 



21k(32k-571-) + 128 



105 (4k3 



n) sm 



4207rK3 



for K < 1 
for K > 1 



(A.2) 



The difference between (A.2| and (A.l) is plotted as a dashed hne in the inset of figure [61 
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